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Abstract 

Feynman diagrams are the best tool we have to study perturbative quantum field 
theory. For this very reason the development of any new technique which allows us 
to compute Feynman integrals is welcome. By the middle of the 80's, Halliday and 
Ricotta suggested the possibility of using negative dimensional integrals to tackle 
the problem. The aim of this work is to revisit the technique as such and check up on 
its possibilities. For this purpose, we take a box diagram integral contributing to the 
photon-photon scattering amplitude in quantum electrodynamics using the negative 
dimensional integration method. The reason for this choice of ours is twofold: Firstly, 
it is a well-studied integral with well-known results, and secondly because it bears 
in its integrand the complexities associated with four massive propagators of the 
intermediate states. 
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Feynman Box Diagrams. 
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1 Introduction. 



Scattering amplitudes, radiative corrections, (3 functions of renormalization 
group, etc., all require the computation of Feynman integrals [1,2], which are 
the more complex to evaluate the more loops in a given diagram one has. 
Since this approach is still the best technique we have to study quantum field 
theory (QFT) perturbatively, solving Feynman integrals becomes basic to do 
any serious study on physical processes of interest involving those quantities. 
Such computations become harder to do not only with increasing number of 
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loops, but also with increasing number of massive particles in the intermediate 
states. This subject becomes even more compeling as we realize that increas- 
ingly higher energies are disposable at particle accelerators and tests for the 
standard electroweak theory are pushed to new limits. It is then clear that 
in order to study radiative corrections there will unavoidably lead us to the 
necessity of dealing with Feynman integrals containing massive intermediate 
vector bosons, and our ability to perform them will play a key role. 

The standard way of solving such integrals starts with the introduction of 
Feynman parameters. Wick rotation and then finally, integration. This method 
is somewhat tedious and sometimes it is not possible to solve exactly the para- 
metric integrals. For this reason physicists developed several other techniques 
to calculate Feynman integrals[3]. A technique known as negative dimensional 
integration method (NDIM) [4-6] has also been considered to tackle the prob- 
lem. 

Our aim in this work is to further check up on NDIM as a useful tool for the 
referred task. For this purpose, we use NDIM to evaluate the box diagram 
contributing to the photon- photon scattering (see fig.l) in quantum electro- 
dynamics (QED). The outhne for our paper is as follows: In section 2 we 
give a brief review of the methodology to be employed, while in section 3 we 
compute the integral proper writing down the two well-known hypergeometric 
series representations for it. In section 4 we consider two particular cases of 
the given integral, namely, integrals with three and two massive propagators 
respectively. Finally, in the last section, we make a few comments about the 
six new results we have for this diagram, which will be the subject to be ad- 
dressed in our shortly forthcoming paper. Also we mention more complicated 
integrals, like the ones arising in one-loop correction to Bhabha scattering in 
QED, and off-shell two-loop three point graphs[7]. 



2 Integration in Negative Dimensions. 



NDIM was introduced by Halliday and Ricotta[4] some years ago. In this 
section we present a brief review of this technique. Basically what one does is 
to perform an analytic continuation 
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Fig. 1. Feynman diagram for photon-photon scattering in the s-channel 

so that one gets a polynomial integral inD < (Qfrom a rather complicated one 
in > 0. We then solve it in D < and go back [5] to D > 0, through another 
analytic continuation. One of the advantages of NDIM is that simultaneously 
we get several hypergeometric series representations for the integral m D > 
0, i.e., we obtain expressions for all the possible regions, physical and non- 
physical alike, of the external momenta. 

We start from the relation[4-6] between a gaussian integral and its counterpart 
in negative dimensions 

/ d^gexp(-Ag^) = {^-^'"' = t^-^j^ / d^Qi^'Y (2) 

where in the last step we have expanded the exponential function in Taylor 
series. Just like in dimensional regularization[9] we take this expression as the 
definition of the negative D— dimensional integral [10]. The middle term is an 
analytic function of D so the integral on the right hand side is also an analytic 

^ The negative dimensional Hnear operator object here can be seen from the view- 
point of positive dimensional fermionic integration [8]. 
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function oi D [11-13]. 

From this equation we get, 

/ d^q{q'y = (-l)%^/%/2,-,r(l + j) 



(3) 



In a similar way, we can solve, e.g.. 



m 



J{i,j,k,l;m)= J d^q (f - 
X (g - k2f' - 
whose counterpart in D > is the integral 



{q — p) —m [q — ki) — m 



(4) 



K{i,j,k,l;m) 



(g^ — m?y [{q — — m?y [{q — kiY — m 
1 



2]k 



[{q - k2Y - m2] 



(5) 



This is one of the integrals that contributes to the photon-photon scattering 
amplitude in QED and it is the one we want to evaluate in our "lab test" 
for NDIM. Of course, since the external photons are real particles, they are 
on-shell, i.e., we consider here that k\ — k2 — {p — k-iY — {p — k2Y — (see 
fig.l). 

So, to begin with, let our "gaussian integral" be 



j d^g exp (-a{q^ - m^) - (3 [{q - pf - 
\q-k2f-m^]) 



m 
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(g — k\Y — 



—UJ 



(6) 



Completing the square, integrating over q and expanding the exponential, we 
get 



ni=0 



ni!n2!n3!n4!n5!n6!n7! 



xr(l +n3-ni-n2 - —) 



(7) 
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where s and t are the Mandelstam variables (see fig.l) and m here is the mass 
of the virtual matter fields. Since we use a multinomial expansion the sum 
index above must satisfy the constraint 

D 

n4 + + riQ + n-i = n3 - Hi - n2 - — 



On the other hand, expanding the exponential of (6), we have 

i,jAi=o r(i + Or(i + j)r{i + k)r{i + 1) 



and comparing the expressions (7) and (8) we obtain a general relation for the 
integral J{i,j,k,l;m) and a system of linear algebraic equations linking the 

sum indices n.j, with five equations and seven variables, i.e., we have twenty- 
one distinct solutions for this system. So, NDIM transfers the problem of 
calculating Feynman integrals to solving systems of linear algebraic equations, 
which is an easier task to perform, of course. 

And for the particular system we are dealing with above, there are six trivial 
solutions, which are of no interest at all and therefore discarded, while five are 
of the hypergeometric type, which will give the known results. The other ten 
solutions — the new results — will be the subject of our shortly forthcoming 
paper. 



3 Hypergeometric Series Representations. 



Solving the system we find five hypergeometric series which can be divided 
into two sets according to its variables, {h}, and {I2, 13, h, h}- The solution 
in the first category is 



X 



t {-i\ni){-j\ni){-k\n2){-l\n2){-a\ni + 



Am?) nM-i-l + f 1^1 + ^2)(| - f + f 1^1 + ^2) 

where we define a — i-\-j-\-k-\-l-\-^ and use the Pochhammer symbol 

( \u\ - ( \ r(a + A;) 
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Substituting i = j — k — I = —1, we get the integral (5) with exponents 
corresponding to the one we want to calculate for box diagrams. Then, the 
first hypergeometric series representation yields (in four dimensions) 



Ji(-l,-l,-l,-l;m) 



6m^ 



1,1, 1,1; 2 



t 



Aw? ' Aw? 



(10) 



This result is exactly what was already obtained by Davydychev[14] using the 
Mellin-Barnes' representation for massive propagators [3]. Note that since we 
are in Euclidean space there is an overall factor i difference when compared 
to Davydychev's result, obtained in Minkowski space (his result has the extra 
factor i). This expression is symmetric in s and t and is nonvanishing for 
s = t = 0. It will allow us, in section 4, to read off the particular cases where 
the integral has three and two propagators respectively, and also the analytic 
continuation to other regions of external momenta. This expression is valid in 
the region of convergence of the series which defines the F3 hypergeometric 
function [15, 16], 

F,(a,a./3,/J,7|x,!,)=E^-^^ ^j-^^^ (11) 



where |x| < 1 and \y\ < 1. In other words, it is valid below the threshold of 
pair production. It is also suitable for studying the non-relativistic limit since 
the Mandelstam variable s must be less than Aw^. We note that s = Aw^ 
defines the point where the process changes its nature, that is, there exists the 
possibility of pair creation and this fact manifests itself in the amplitude as a 
branch point in the Feynman integral [2, 13]. 

The hypergeometric function in (10) can be expressed in terms of its dou- 
ble integral representation [15], and from there one can arrive at the standard 
result expressed in terms of a rather cumbersome sum of logarithms and dilog- 
arithms[14,17]. 

The next set of solutions, also obtained by Davydychev[14], has variables 
Aw^ / s and Arr?/t, the inverse of the ones in the first solution. In the following 
we write down these four solutions: 



/2 = 



2n^/^{-ty{-sy{-w^)D/^+'+''{-i\j){-k\l) 



{l+i-j + k-l\§+j + l) 



E 

ni, 712=0 



Am 



2\ni 



ni\n2\ 



Am 



2\ "2 



X 







"2) 


^ i+i-i+k-j 


III + 11-2 


) 




(l + i- j|ni)(l + A;-/ 


n2)(l + i + A; + f 1 


ni + ^2) 



X 



(l + |(i-j+A;-0|ni + n2), 



(12) 
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27r^/2(-i)'=(-sy(-m2)f+*+'(-i|j)(-/|A;) 



(l + i + j-k-l\§+j + k) 



E 



4m 



X 



"2 



(-i|nO(-A;|n2) ( 



ni,n2- 

1+i-j-k+l I 
2 



=0 niln2l 



t J (l + i-j|ni)(l-fc + %2)(l + i + i + f 1^1+712) 



X (l + i(i - j - A; + /) 



ni + 77,2 



(13) 



'4m2r^ (-7|ni)(-A;|n2) ( i-^+W | ^ 



and 



X 



X 



V t J (l-i+j|ni)(l-A; + /|n2)(l+j + / + f|ni+n2) 

{l + l{-i+j-k + l)\ni + n2), (14) 



^ _ 27r-P/^(-t)'(-5)^(-m^)f +i+fe(-j|i)(-fc|/) ^ 1 Am^"^ 



X 



t J {l-i+j\ni){l + k-l\n2){l+j + k + ^D\ni + n2) 
X (l + l{-i+j + k-l)\ni + n2) (15) 

Prom these we construct the second hypergeometric series representation for 
the Feynman integral as just the linear combination 72 + -^3 + -^4 + -^5, i-e., 



M,j.k,,,,„)- ^L^^ ^^^^^^^^ (16) 



XF2 (l + a + 6 — 7,a,6;l + a — c, 1 + 6 — (i 



4?TT,^ 4m^' 



t 



where the set {a, b, c, d} takes the values {—i, —k, —j, —I}, {—i, —I, —j, —k}, 
{—j, —k, —i, —I}, {—j, —I, —i, —k} and F2 is another hypergeometric func- 
tion[15,16]. 

Note that J2 is the analytic continuation of Ji, see for example Erdelyi et al 
[15]. While Ji is valid in the non-relativistic case and J2 in the relativistic one, 
they do not cover all the possible regions of external momenta. 

An important point here that one must be aware of is that even though singu- 
larities might appear in isolated terms of the RHS, for the special case when 
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i — j — k — I = —1, the above equation (16) cannot be singular since the 
corresponding analytic continuation formula (10) is not. 

To overcome this difficulty let us introduce small corrections in the parame- 
ters P and P' of the hypergeometric function F3{a,a', P, P';^\x,y)[18]. This 
recourse corresponds to correcting the exponents of propagators [19]. In our 
case we take 

P'^l + 5'. 



Next we expand all the (gamma) factors which contain 6 and 6' and the F2 
functions in Taylor series around 6 = and 5' = 0. In the end we take the 
limit of vanishing 5 and 5'. This result is valid, like in the first case, within 
the region of convergence of the series which defines the F2 function [15, 16], 

F2{a,P,P;^,-f\zi,Z2) = 2^ —jj , , , ^1 + ^2 < 1 



Following these steps we arrive at Davydychev's second expression for the 
Feynman integral [14]. Note that now the Mandelstam's variables s and t can 
never be zero. Moreover, like the case of Ji, he has shown that the resulting 
expression for J2 can be converted into the standard result in terms of sum of 
logarithms and dilogarithms through the use of a double integral representa- 
tion for F2[15]. 

We see that with NDIM we immediately get two hypergeometric functions 
which are related by analytic continuation. But as we have mentioned before, 
there are still ten other solutions for the system, so that we expect that these 
new results will generate solutions for the integral that are also related by an- 
alytic continuation. The theory of generalized hypergeometric functions[15,20] 
tells us that there are three analytic continuations with no restriction in the 
parameters, namely: 

Fs{...\x,y) ^ F2{...\l/x,l/y). 



This is the very relation we saw above, connecting the two hypergeometric 
functions. The other two, 

Fs{...\x,y)^H2{...\l/x, -y) 
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and 



H2{...\x,y)^F^{...\l/x, -y) 



must appear in our new results whether in the form of direct or indirect 
analytic continuation. So these new results for the Feynman integral (5) will 
cover all the possible physical regions of external momenta. 



4 Particular Cases. 



If we put in our integral any exponent equal to zero we get an integral with 
three massive propagators. NDIM must be consistent for any value of the 
exponents, so, taking for example i = the result (9) gives 



Ji(0,j,A;,Z;m)=C« (17) 



/ — fc, — —J — k — I — , 

3-C^2 I -j-k-l l-j-k-l 

\ 2 ' 2 



t 



Arn? 

where 3F2 is the generalized hypergeometric function of one variable[15] and 

fr(-j-fc-z-f) 



r7« = (—\'" 2v/^(-2^')^'^'^'^^r(-^ 
V 2 y r(^i^)r(i^i^) 



The same result arises if one puts j = 0. For A; = we need simply to replace 
s f to get 



Ji(^,j,0,/;m)=C(2) (IJ 

s \ 



I, I J I 

i— 

2 ' 2 



where 



2 ; r(=^)r(i=^) 



Finally, if we take / = we get the same result as above in virtue of the 
symmetry {k <-> I) in equation (9) . 
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Prom these equations above we can also calculate integrals with two massive 
propagators, i.e., vacuum polarization-like ones. In (17) taking for example 
j = 0,we get 

^ 2 ' 2 




where 



V 2 ; r_2m2)-*=-'-f r(^)r(i^) 



In an analogous way, if we take / = in (18), we have 



Ji(z,j,0,0;m) = C(^)3F2( dl-j ' ^ 



2 ' 2 



D 

4m^ 



where 

D 



^(4) ^ (zlY^' 2Vir(-z-j-f) 
V 2 ; (_2m2)-i-i-f r(^)r(i^) 



For the particular cases where the exponents are equal to minus one the results 
above are respectively given by. 



Ji(-l,-l,0,-l;m) 



-TT 



2m2 



3-C^2 



1, 1, 1, 

2 ^ 

^' 2 



4m^ 



(19) 



which can be written in terms of elementary functions[3] . Here we have already 
taken the limit D — A, since the result is finite in this limit, and 



Ji(0,0,-l,-l;m) =7r2(m2) 



2 - 



(20) 



where now, of course, there is a well-known single pole in the limit = 4. 

If one is interested in other regions of external momenta one has to take the 
analytic continuation [21] of these results. We do not need to start from the 
very beginning, eq.(9), to construct them in other regions, because we know 
from the theory of complex variables that any functional property of one power 
series is shared by all the others [2 2]. 
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5 Conclusion. 



Using NDIM we have evaluated a massive box diagram integral, namely, a 
Feynman integral bearing four massive propagators. This integral is the one 
appearing in the photon-photon scattering in QED and the two well-known re- 
sults, expressed in terms of hypergeometric functions, have been easily found. 
So, the computation of such an integral, done as a "lab test" for NDIM, has re- 
vealed to us a powerful technique, which transfers the intricacies of performing 
Feynman integrals in positive dimensions to that of solving a system of linear 
algebraic equations in negative dimensions, a far simpler task to perform than, 
e.g., solving parametric integrals. More than that, surprisingly, the technique 
not only reproduces the standard results, but gives simultaneously, solutions 
covering other regions of the external momenta. We are studying carefully such 
solutions and these are going to be the subject of our shortly forthcoming pa- 
per. Also more complicated box diagram integrals, such as arising from the 
one-loop correction to the Bhabha-scattering in QED and an off-shell two-loop 
triangle diagram integral have already been checked by us. 
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